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\ Abstract 

The transverse electromagnetic waves propagating in a relativistic two-fluid plasma influ- 
^ ■ enced by the gravitational field of the Reissner-Nordstrom-de Sitter black hole has been 
investigated exploiting 3-1-1 split of spacetime. Reformulating the two-fiuid equations, 
^ ■ the set of simultaneous linear equations for the perturbations have been derived. Using 
a local approximation, the one-dimensional radial propagation of Alfven and high fre- 
[ quency electromagnetic waves are investigated. The dispersion relation for these waves is 
^ obtained and solved numerically for the wave number. 
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1 Introduction 



Black holes belong to the most fascinating objects predicted by Einstein's theory of grav- 
itation and are still mysterious [1]. Physicists are grappling the theory of black holes, 
while astronomers are searching for real- life examples of black holes in the universe [2]. 
However, the theory of general relativity and its application to the plasma close to the 
black hole horizon have remained esoteric, and little concrete astrophysical impact has 
been felt. Within 3Rs (3 Schwarzschild radii) it is possible to have plasma [31 IH El E]- 
The plasma in the black hole environment may act as a fluid and black holes greatly 
affect the surrounding plasma medium (which is highly magnetized) with their enormous 
gravitational fields. Hence plasma physics in the vicinity of a black hole has become a 
subject of great interest in astrophysics. In the immediate neighborhood of a black hole 
general relativity applies. It is therefore of interest to investigate plasma physics problems 
in the context of general relativity. 

In recent years, considerable attention has been concentrated on the study of black 
holes in de Sitter (dS) spaces. The motivation behind it is based on two aspects: first, 
several different types of astrophysical observations indicate that our universe is in a 
phase of accelerating expansion 0ll[9llIOllIIlliail3l[Tl[T5l[T6l[I7l[l8l[^ 
[22| [231 [211 [251 [261 [23 [281 [29], associated with which is a positive cosmological constant. 
Our universe therefore might approach a dS phase in the far future [301 ISH [321 [331 [Ml 
[351 [36] . Secondly, like the AdS/CFT correspondence, an interesting proposal, the so- 
called dS/CFT correspondence, has been suggested that there is a dual relation between 
quantum gravity on a dS space and Euclidean conformal field theory (CFT) on a boundary 
of dS space [371 [2H1 [391 [iQl [ID [12] . In view of these reasons, the study of transverse wave 
propagation in relativistic two-fluid plasma in the vicinity of a black hole in dS space is 



A covariant formulation of the theory based on the fluid equations of general relativity 
and Maxwell's equations in curved spacetime have so far proved unproductive because 
of the curvature of four- dimensional spacetime in the region surrounding a black hole. 
Thorne, Price, and Macdonald (TPM) [l3l [lH [13 [16] developed a method of 3+1 formu- 
lation of general relativity in which the event horizon of black holes was replaced with a 
membrane endowed with electric charge, electrical conductivity, and flnite temperature 
and entropy so that the physics outside the event horizon turns out to be very much sim- 
pler than it would be using the standard covariant approach of general relativity. Buzzi 
et al. [171 [IH] employed the TPM formulation to develope a general relativistic version of 
two-fluid plasma and investigated the nature of plasma in the vicinity of Schwarzschild 
black hole. Exploiting the formalism of Buzzi et al., we investigate the transverse wave 
propagation in the two-fluid plasma around the Reissner-Nordstrom-de Sitter (RNdS) 
black hole which is the Schwarzschild black hole generalized with a charge parameter and 
a positive cosmological constant, described by the metric 



interesting. 



ds'^ = -f{r)dt^ + 



+ r\d9^ + sm^9d^^), 



(1) 



fir) 



where the metric function /(r) is given by 



fir) = 1 



2M g2 




(2) 



2 



M being the mass, and q the total charge (electric plus magnetic) in the background 
of the static de Sitter space. The positive cosmological constant A is usually written as 
A = 3/£^ with £ the cosmological radius. The spacetime causal structure depends strongly 
on the singularities of the metric given by the zeros of f{r). Depending on the parameters 
M, g, and A, the function /(r) may have three, two, or even no real positive zeros. For 
the RNdS black hole case we are interested in which /(r) has three real, positive roots 
(r_ < r+ < Tc), and a real negative root r„ = — (rc + r-|. + r„). The horizons at r_, r_|_, and 
Tc, are called the inner (Cauchy), outer (event), and cosmological horizons, respectively. 
The black hole parameters M, g, and £ are related to the roots by 



r'\ + r1 + r^Tc — r_rn 
(r+ + rc){r+rc - r_r„) 



2fM, 



For the zeros of f{r), we have 

- fr^ + 2Mfr - fq^ = 0. 
Positions of the black hole horizons are given by 



(3) 
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The cosmological horizon is located at 

/ 



sm 



1 



sm 



3V3M 



1 + V 



1 + 



1 + c SVSM 



CSC-' 



sm 



J- -h ^2 



- 1 



(7) 



Since 27^ < 1, therefore 3^3^^/ VI + ^ < 1, and we obtain 
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It gives the the RN black hole horizons for i ^ oo and the Schwarzschild-de Sitter black 
hole horizon for q = 0. The metric ([T]) represents an interesting asymptotically de Sitter 



extreme RN black hole 



> 



i+C 



it 



(or, "cold" black hole) for q^ = j^M , while for q 
does not represent any black hole but an unphysical naked singularity. 

The black hole horizon r_|_ and cosmological horizon are not in thermal equilibrium 
because the time periods in the Euclidean section required to avoid a conical singularity at 
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both do not match in generaL However, there are two famihes of RNdS solutions with the 
same black hole and cosmological horizon temperatures, representing the "lukewarm" black 
hole and charged Nariai black hole [IHl ED] • The lukewarm black holes are characterized 
by = ^'^1 which is the condition of stable endpoints of the evaporation process. Its 
three horizons do not coincide and cosmological horizon and black hole horizon satisfy 

(SUES] 

M^ = q^ = -i^ E = r^ + r,. (10) 
(r+ + r^Y 

The black hole is colder than the lukewarm solution for > and it will then absorb 
radiation from the cosmological horizon. As increases relative to the inner and outer 
horizons eventually coincide. The resulting hole then becomes extremal (or "cold") RNdS 
black hole with the mass, charge, horizon radius and cosmological constant satisfying the 
following relations |52j : 

M ,^11 



r_|_(rc + r+)^ r^r'^lrc + 2r^) i'^ + 2r+rc + 3r^ 

As M increases relative to |g| with M > 0, the outer black hole and cosmological 
horizons come closer together. When these horizons coincide at Vh, where f{rh) = f'{fh) = 
0: 

- - > ^ fW) ' (^^' 

it then becomes the charged Nariai solution, which is the largest charged asymptotically 
dS black hole. For a given q, it has the maximal mass Mn- The charged Nariai solution 
and lukewarm solution join together for the critical value q^ = = i'^/16. If M^m is 
the mass of extremal RNdS black hole, the range for mass parameter of a black hole in 
dS space is M^m < < Mjy, referred to as the "undermassive"case. The metric ([1]) 
describes a naked singularity, if this limit is exceeded. 

If the three horizons coincide, i.e. = r_|_ = r_, the RNdS black hole becomes an 
"ultracold" black hole with horizon at r^cd, where fivucd) = f'i^ucd) = f'ijucd) = 0. The 
mass, charge and cosmological constant are related according to 

M = h'ucd, q' = \rl^,, f = Qrl,. (13) 

This configuration simultaneously maximizes the values of M and for any given positive 
value of the cosmological constant A. 

The organization of this paper is as follows: In section [2], we summarize the 3+1 
formulation of general relativity and Maxwell's equations in curved spacetime. In section 
[21 we describe the horizon governing equations of two-fluid plasma in the RNdS black hole 
spacetime. In section IH we consider the wave propagation in the radial z direction, and 
linearize the equations for wave propagation in section [5] by giving a small perturbation 
to fields and fluid parameters. In section [6|, we derive the dispersion relation for the 
transverse waves and discuss the local approximation used to obtain numerical solution 
of it. In section [TJ we discuss our numerical solution modes and present the results in 
section [HI Finally, we give our concluding remarks in section [3 Throughout the paper we 
use natural units: G = c = = 1- 



4 



2 3+1 Formulation in Spacetime 



An absolute three-dimensional space defined by the hypersurfaces of constant universal 
time t is described by the metric 

ds^ = gijdx^dx^ = JJ^^^'^ + r'^idO'^ + sin^ Odip"^). (14) 

The indices i,j range over 1,2,3 and refer to coordinates in absolute space. The fiducial 
observers (FIDOs), that is, the observers at rest with respect to this absolute space, 
measure their proper time r using clocks that they carry with them and make local 
measurements of physical quantities. Then all their measured quantities are defined as 
FIDO locally measured quantities and all rates measured by them are measured using 
FIDO proper time. The FIDOs use a local Cartesian coordinate system with unit basis 
vectors tangent to the coordinate lines: 

d Id Id 



For a spacetime viewpoint rather than a 3+1 split of spacetime, the set of orthonormal 
vectors also includes the basis vector for the time coordinate given by 

did , ^ 

dr a at 

where a is the lapse function (or redshift factor) defined by 

, , dr / 2M g2 ^2\| 

The gravitational acceleration felt by a FIDO is given by 

1 f M r\ 



-Vina = - 



while the rate of change of any scalar physical quantity or any three-dimensional vector 
or tensor, as measured by a FIDO, is defined by the derivative 

D ( I d \ 
V being the velocity of a fluid as measured locally by a FIDO. 



3 Two-fluid Plasma Governing Equations 

We consider two-component plasma such as an electron-positron or electron-ion. In the 
TPM notation, the continuity equation for each of the fluid species is given by 

d 

Olilsris) + V ■ {a-isUs^s) = 0, (20) 
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where s is 1 for electrons and 2 for positrons (or ions). For a perfect relativistic fluid 
of species s in three-dimensions, the energy density e^, the momentum density S^, and 
stress-energy tensor W^'' are, respectively, given by 

6, = 7.'(es + Ps^l), S, = 7.'(^s + Ps)v„ = 7,^(5, + PsHv'^ + P,g^■^ (21) 

where is the fluid velocity, is the number density, Pg is the pressure, and Eg is the 
total energy density defined by 

Es = msUs + Ps/ilg - 1). (22) 

The gas constant 7g is 4/3 for T — > cxd and 5/3 for T — > 0. 

When the two-fluid plasma couples to the electromagnetic fields. Maxwell's equations 
take the following form: 

V-B = 0, (23) 
V-E = 47ia, (24) 

— = -Vx(aE), (25) 

— = V X (aB) - 47raJ, (26) 
where the charge and current densities are respectively defined by 

s s 

The fluid quantities in (l2TI) take the following form in the electromagnetic field: 

e, = ^{E' + B'), S. = -^ExB, 

Wf = — (E^ + B^)g^^ - —{E^E^ + B^B''). (28) 
Stt 47r 

Energy and momentum conservation equations are respectively expressed by ^31 HH US], 

-|:e. = -V-S, + 2a-S„ (29) 
a at 

-|:S, = e,a--V-(aW.), (30) 
aot a 

which, using (!22l) - (!26l) . can be written for each species s in the form 

+7,g,n,E ■ v, + 2^'^^{Es + Ps)a ■ v, = 0, (31) 
7,^(es + Ps) I + V, ■ V 1 V, + VPs - IsQsnsiE + V, X B) 



a dt 

a dt 



+v, (-fsqsnsE ■ Vs + -^Ps] + Isi^s + ^.)[v.(v, ■ a) - a] = 0. (32) 



Although these equations are vahd in a FIDO frame, they reduce to the corresponding 
special relativistic case for a = 1 [33], which are valid in a frame in which both fluids are 
at rest. The transformation from the FIDO frame to the comoving (fluid) frame involves 
a boost velocity, which is simply the freefall velocity onto the black hole, given by 

vs = {l-a')l (33) 

Then the relativistic Lorentz factor 7boost = (1 — ^s)^^ = 1/a. 

The ion temperature profile is closely adiabatic and it approaches 10^^ K near the 
horizon [54j. Far from the (event) horizon electron (positron) temperatures are essentially 
equal to the ion temperatures, but closer to the horizon the electrons are progressively 
cooled to about 10*^ — 10^ K by mechanisms like multiple Compton scattering and syn- 
chrotron radiation (i.e. radiation which occurs when charged particles are accelerated in 
a curved path or orbit). Using the conservation of entropy the equation of state can be 
expressed by 

^(P>lg)=0. (34) 
Dt 

where D/ Dt = {l/a)d/dt + v,, • V. The full equation of state for a relativistic fluid, as 
measured in the fluid's rest frame, is given by [551 EH] 



msUs + VflsHs 



Ps iH^^^'iimsns/Ps) 
iH^\imsns/Ps) 



(35) 



where the H'2\x) are Hankel functions. 

The Rindler coordinate system, in which space is locally Cartesian, provides a good 
approximation to the Reissner-Nordstrom-de Sitter geometry near the event horizon. The 
essential features of the horizon and the 3+1 split of spacetime are preserved without the 
complication of explicitly curved spatial three-geometries. In Rindler coordinates the 
Reissner-Nordstrom-de Sitter metric is approximated by 

ds^ = -a^dt^ + dx^ + dy^ + dz^, (36) 

where 

X = r+{e - ^ , y = r+ip, z = 2r+a. (37) 

The standard lapse function simplifies in Rindler coordinates to a = z/2r_|_, where r+ is 
the event horizon of the black hole. 



4 Wave Propagation in Radial Direction 

We consider one-dimensional wave propagation in the radial z direction. Introducing the 
complex variables 

Vsz{z, t) = Us{z, t), Vs{z, t) = Vsx{z, t) + Wsy{z, t) , 

B{z,t) = B^{z,t) + iBy{z,t), E{z,t) = E^{z,t) + iEy{z,t), (38) 
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we obtain 



and defining 

B, = ^{B + B*), By = j^{B-B*), (39) 

Vs,By - VsyB^ = ^{v,B* - vlB), (40) 

and similarly for the transverse components of E, where the * denotes the complex con- 
jugate. The gravitational acceleration becomes 

d 1 da 

The equation of continuity ( 120|) and Poisson's equation (HM are respectively written 

as 

d d 

g^h'^^'^) + g^iajsnsUs) = 0, (42) 

and 

= 47r(gini7i + 52^272)- (43) 
From the and components of (|25|) and fl26|) . one could derive 

if - - 

dE f d \ 

i— = -a{ — -a\B- i47rea (72^2^72 - 7i^i^i)- (45) 
Differentiating (145!) with respect to t and using (jH]), we obtain 

"V + '^^I^^-^^+i^jj j^ = 4-ea^(-272.2-n,7i.i). (46) 

From the and &y components of fl32l) . the transverse component of the momentum 
conservation equation is obtained in the form 

Ps^^ = qsUslsiE - iVsB^ + lUsB) - u^Vspsa - Vs ( QsUs'^s^ ■ ^.s + --gf 1 ' (47) 
where 

and ps is the total energy density defined by 

Ps = -fUss + Ps) = 7'(^sn. + FgP,), (48) 

with Fg = 7g/(7g - 1). 
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5 Linearization 



We use perturbation method to linearize the equations derived in the preceding section. 
We choose an apphed magnetic field to lie along the radial direction and introduce the 
quantities 

Us{z,t) = uos{z) + 6us{z,t), Vs{z,t) = 6vs{z,t), 
nsiz,t) = nosiz) + Snsiz,t), Ps{z,t) = Pq,{z) + 5Ps{z,t), 
psiz,t) = posiz) + 6psiz,t), E{z,t) = 6E{z,t), 
B,{z,t) = Boiz)+6B,{z,t), B{z,t) = 6B{z,t). 



(49) 



The relativistic Lorentz factor is also linearized such that 



7. 



70s + S'js, where 7os = (1 



u 



Os 



(50) 



The unperturbed radial velocity near the event horizon for each species as measured by 
a FIDO along is assumed to be the freefall velocity so that 



Uos{z) = vs{z) = [1 - a'^iz)] 
which, with ((31) and ( fT7|) . gives 



vs{z) 



X, 



(51) 



(52) 



where 



X 



1 + 



2M 



1 + — — - — 

2M\i^ rl, 



(53) 



For the plasma falling at the black hole event horizon, vq^ = 1 and at cosmological horizon, 
t>g = (r_|_/rc) < 1, while for the plasma inside the cosmological horizon and falling at the 
black hole horizon, we have 1 < x < {'^+/^cY^'^- It follows, from the continuity equation 
fH2]) . that r"^ d'-^Qs'n'OsUQs = const. = r^Q;_|.7+'n,+u+, where the values with a subscript 
"+"are the limiting values at the event horizon. The freefall velocity at the horizon 
becomes unity so that = 1. Since Uqs = "^^fr, 7os = and hence a7os = «+7+ = 1- 
Also, because = x{'^+I^Y^'^ i the number density for each species can be written as 
follows: 

1 



n^sKz) = —n+sVfi[z. 
X 



(54) 



The equation of state (IMl) and (l5ll) lead to write the unperturbed pressure, in terms of 
the freefall velocity, as follows: 



Pos{z) 



1 



370 



X 



{z). 



Since P{ 



Os 



kBUosTos, then with fc^ = 1, the temperature profile is 

1 „ 3(7g-l), 



Tosiz) 



X 



{z). 



(55) 



(56) 
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The unperturbed magnetic field is purely in the radial direction. It does not experience 
effects of spatial curvature. From the flux conservation V-Bq = 0, it follows that r'^Bo{r) = 
const., from which one obtains the unperturbed magnetic field, in terms of the freefall 
velocity as 

Bo{z) = -,B^vUz), (57) 
X 

where vs{z) = [1 — 0^(2;)]^/^. Since 

^„ , ^,1 

(5^ 

az ZT+ vs 

we have 





dvs 


a 1 






dz 


2r+ vs ' 




duos _ 


a 1 


dBo 


4a Bq 


dz 


2r+ vs ' 


dz 


2r+ vl ' 


driQs _ 


3a nos 


dPos _ 


3a 7gPos 


dz 


2r+ vl ' 


dz 


2r+ vl 



(59) 

When the linearized variables from (H9l) and ( |50l) are substituted into the continuity 
equation and products of perturbation terms are neglected, it results 

f d d da 2 duQs\ r ( d da\ , , 

[-di + + + j ^[^d-z^Tz) 



d d da / 1 dnQs ^ du^s 
The conservation of entropy, when linearized, gives 



5m, = 0. (60) 



= ^5n., (61) 



and hence from the total energy density, fj48l) . it is possible to write 



= — ( 1 + '^^ 1 <5^s + '2Uqs11sP0s8Us, (62) 



where pos = 7os(^s^Os + rg-Pos)- When the transverse part of the momentum conservation 
equation is linearized, differentiated with respect to t, and substituted from (jH]), it gives 

/ d d da iaqsjosnosBo\ d5vs 



'dz dt dz J dt 

aqslosUos ( d d 9a\ 

aUosTT- + 777 + Uos^ ] oE = 0. (63) 



pos \ dz dt dz 
Poisson's equation ( H3l) and ( H6l) are linearized to obtain, respectively, 

^ = 47re(no27o2 - '^oiToi) + 47re(7o25n2 - loiSui) 

+47re(no2Mo27o2'^W2 - ?^oiMoi7oi'^mi), (64) 

2 9^ 9a; 9 (9^ f da\\ ( d5v2 ddvA 

" ^ + 7^-7^+ ( — ) |5^ = 47rea(no27o2^-rioi7oi^)- (65) 



da d d'^ f da\\ ( d5v2 ddv] 

'Tzd-z'W^[-d-z) j = 1,^°^^°^^ - 
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6 Dispersion Relation 



We investigate the information which may be obtained from the hnearized equations 
derived above by restricting our consideration to effects on a local scale for which the 
distance from the horizon does not vary significantly. We use a local (or mean-field) 
approximation for the lapse function a{z) and hence for the equilibrium fields and fluid 
quantities. Since the plasma is situated relatively close to the horizon, a relatively small 
change in distance z will make a significant difference to the magnitude of a. It is thus 
important to choose a sufficiently small range in z for which a{z) does not vary much. 
We consider thin layers in the direction, each layer with its own ao, where ao is some 
mean value of a within a particular layer. Then a more complete picture can be built up 
by considering a large number of layers within a chosen range of ao values. 

The local approximation imposes a restriction on the wavelength and hence on the 
wave number k. It is assumed that the wavelength is small compared with the range 
over which the equilibrium quantities change significantly. Hence, the wavelength must 
be smaller in magnitude than the scale of the gradient of the lapse function a, i.e.. 



A < 

or, equivalent ly, 

2n 




1r, 



1 

2^ 




2e 

1 + CM2 



5.896 X lO^cm, 



k > 



2r, 



2e 




1 + CM2 



1.067 X 10"W"\ 



with 1 < 



1 + C 



< 



g2 



for a black hole of mass ~ IMq. The value for 2^/(1 + (^) = 1 corresponds to the RN 
black hole, while the value for 2^/(1 + (^) = M^/g^ corresponds to the extreme RNdS 
black hole. 

Since the hydrodynamical approach used in the work is essentially a bulk, fluid ap- 
proach, the microscopic behavior of the two-fluid plasma is treated in a somewhat approx- 
imate manner via the equation of state. Then the results are really only strictly valid in 
the long wavelength limit. However, the restriction, imposed by the local approximation, 
on the wavelength is not too severe and permits the consideration of intermediate to long 
wavelengths so that the small k limit is still valid. 

The unperturbed fields and fluid quantities are not assumed to be constant with respect 
to a{z). Then, the derivatives of these quantities can be evaluated at each layer for a 
given ao- In the local approximation, a ~ ao is valid within a particular layer. Hence, 
the unperturbed fields, fluid quantities and their derivatives take on their corresponding 
"mean- field" values for a given ao. The coefficients in f l60p . (1631) . and (1641) become constants 
within each layer, evaluated at each fixed mean- field value, a = ao. So, it is possible to 
Fourier transform the equations with respect to z, using plane- wave-type solutions for the 
perturbations of the form ~ ^i{kz-ujt) each ao layer. 

When Fourier transformed, (jH]), fHSl) . fl63|) . and fl65l) turn out to be 



k —\5E + —5B = Q, 

. 2r+ao / ao 



'^6E=(k 
ao V 



2r+ao^ 



5B + 47re(7o2no25t;2 - 7oi"'Oi'5^i) 



(66) 
(67) 
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and 



/ , , iuos , aoqs'yosnosBo] . 

LU aokuos -uj + \ dvs 

V 2r+ pos J 

. Qslosnos ( , \ . p „ 

-lao dQkuQs -uj - - — = 0, 

Pos \ 2r+ ; 



5E 



(69) 



aok{aok - i3/(2r+)) - uj^ - l/(2r+)2 ' 

The dispersion relation for the transverse electromagnetic wave modes may be put in the 
form 



an 



2r^ 



1 



(2r+)2 
+ 



M02-^=l 



1 (Woi-f^T ~bJ - ao^^cl) (%2-^T - + "0'^c2) 



(70) 



for either electron-positron or electron- ion plasma, where K± = a^k ± i/2r_| 

e^7os'n.os-Bo/pos and cOps = ^j47i^^^y^~n^Jpo^. Like the plasma frequency cUpg, the cyclotron 
frequency Ucs is frame independent. The fluid quantities are measured in the fluid frame, 
but the field Bq is measured in the FIDO frame. Hence, the factors of 70s do not cancel 
out explicitly. The transformation Bq — > 70^-60 boosts the fluid frame for either fluid 
and thereby cancels the 70s factors. The "+"and "—"denote the left (L) and right (i?) 
modes, respectively. The dispersion relation for the L mode is obtained by taking the 
complex conjugate of the dispersion relation for the R mode. The two modes have the 
same dispersion relation in the special relativistic case. 





Figure 1: Left: Real part of Alfven damped mode for the electron-positron plasma. Right: 
Imaginary part of the damped mode. 



7 Numerical Solution of Modes 

The dispersion relations ( 1701) are complicated enough to make any attempt at an analytical 
solution even in the simplest cases for the electron-positron plasma where both species 
are assumed to have the same equilibrium parameters. As the analytical solution is 
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Figure 2: Top: Real part of the complex conjugate pair of Alfven damping and growth 
modes for the electron-positron plasma. Bottom: Left: Imaginary part of the damped 
mode. Right: Imaginary part of the growth mode. 



cumbersome and unprofitable, we solve the dispersion relation numerically and for this 
purpose we put it in the form of a matrix equation as follows: 

{A - kI)X = 0, (71) 

where the eigenvalue is chosen to be the wave number k, the eigenvector X is given by 
the relevant set of perturbations, and I is the identity matrix. 

We need to write the perturbation equations in an appropriate form. The equations 
are then written in terms of the following set of dimensionless variables: 

^ cu ~ kc 1 
oj — , k — — , /c+ = , 

SUs . SVs SUs 
dUs — , Vs — , OUs = , 

uos Uf)s rios 

-Do -Do -Do 

For an electron-positron plasma, Upi = u}p2 and uid = u!c2', hence, a;* is defined as 

{Uc Alfven modes, 

{2u}l + u}'^)^ high frequency modes. 
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Figure 3: Top: Real part of the complex conjugate pair of Alfven damping and growth 
modes for the electron-ion plasma. Bottom: Left: Imaginary part of the damped mode. 
Right: Imaginary part of the growth mode. 



where ujp — y/uJ^iUJ^ and lUc = ■\J^c\^c2- However, for the case of an electron- ion plasma, 
the plasma frequency and the cyclotron frequency are different for each fluid. The choice 
of Lo^ is then a more complicated matter. For simplicity, we assume that 



a;* 



+ a;^2) ^ Alfven modes, 
^Ix + ^1-2)^ high frequency modes. 



(74) 



where = {2Upg + u'^g). In the zero gravity hmit, these values for reduce to the 
special relativistic cutoffs for electron-positron plasma. The special relativistic cutoffs 
are determined by the dispersion relation in view of the fact that the solutions to the 
dispersion relation arc physical (i.e., Re(A;) > 0) only for certain frequency regimes. As 
the dispersion relation can not be handled analytically, it is difficult to determine the 
cutoffs in the case including gravity. Other similar combinations for a;* make only a 
difference of a scale factor to the form of the results as a;* is really only a scale factor. 
We write the dimensionless eigenvector for the transverse set of equations in the form 



transverse 



6B 
SE 



(75) 
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Figure 4: Left: Real part of Alfven growth mode for the electron-ion plasma. Right: 
Imaginary part of the growth mode. 



Using (I72|) . we write flBB]) - flB5]) in the dimensionless form as follows: 

~ ik ~ 

k6E = -iu6B + —6E, (76) 
k6B = uoi^^6vi - uo2^^Sv2 + —SB + iu6E, (77) 

and 




.SB-ii^] ^^6E. (78) 
Mosto'* ao y V e / uqs^^ \ e J uqsu;* 

These are the equations in the required form to be used as input to flTT]) . 

8 Results 

We have carried out the numerical analysis by using the well known MATLAB. We have 
taken the values = 0.5M^, x = 0.9, and ( = 1.8. We have considered both the 
electron-positron plasma and the electron-ion plasma. The limiting horizon values for the 
electron-positron plasma are taken to be 

n+s = 10^^ cm-^ T+, = 10^° K, 5+ = 3 x 10^ G. (79) 

For the electron-ion plasma, the ions are essentially nonrelativistic, and the limiting hori- 
zon values are chosen to be 

= lO^s cm-3, T+i = 10^° K, n+2 = 10^^ cm'^, T+2 = 10^^ K. (80) 

The equilibrium magnetic field has the same value as it has for the electron-positron 
case. The limiting horizon temperature for each species has been taken as derived by 
Colpi et al. [51] from studies of two temperature models of spherical accretion onto black 
holes. For the limiting horizon densities and the limiting horizon field, simply arbitrarily 
those values have been chosen which are not inconsistent with the current ideas. The gas 
constant and the mass of the black hole have been chosen as follows: 

4 

7ff = - and M = 5Mq. (81) 
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Figure 5: Top: Real part of the growth mode and of the damping and growth mode for 
the high frequency transverse wave in the case of the electron-positron plasma. Bottom: 
Left: Imaginary part of the growth mode. Right: Imaginary part of the damping and 
growth mode. 

8.1 Alfven Modes 

8.1.1 Electron-positron Plasma 

For electron-positron plasma four modes exist, two of which are growth and other two 
are damped. The first mode, shown in Fig. 1, is damped. The second mode, not shown 
here, is equal to the first mode but opposite in sign and is growth. The third and fourth 
modes, shown in Fig. 2, are complex conjugate pairs and show damping and growth, 
respectively. These two modes are similar to the corresponding Alfven modes of Buzzi 
et al. [U] but the damping and growth rates are larger by several orders of magnitude 
compared with those modes. These four modes coalesce with two modes of Buzzi et al. 
|17] for the Schwarzschild case and with a single mode for the ultra relativistic electron- 
positron plasma [23] • Since we are using the convention e'*^^ = ei[Re(fc)+iim(fc)]z^ ^^^q damping 
corresponds to Im(k) > and growth to Im(k) < 0. 

8.1.2 Electron-ion Plasma 

In the case of the electron-ion plasma there exist four modes, two of which, shown in Fig. 
3, are complex conjugate pair and are significantly damped and growing respectively. 
The other two modes are also damping and growth. The third mode, shown in Fig. 4, 
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Figure 6: Top: Real part of the damped mode and of the damping and growth mode for 
the high frequency transverse wave in the case of the electron-positron plasma. Bottom: 
Left: Imaginary part of the damped mode. Right: Imaginary part of the damping and 
growth mode. 

is a growth mode. The fourth mode is similar to this mode but opposite in sign and is 
damped. The last two modes display only marginal damping and growth, respectively, 
and are equivalent to the two modes discussed in Fig. 2 for electron-positron plasma. For 
the first two modes, shown in Fig. 3, the differences in the magnitudes of the Ud and 
uJc2 apparently lead to take the frequencies from their negative (and therefore unphysical) 
values for the electron-positron case to positive physical values for the electron-ion case. 
These changes are because of the difference in mass and density factors as between the 
positrons and ions. These four modes coalesce with the modes of Buzzi et al. |17] for 
the Schwarzschild case. It is evident that the growth and damping rates are independent 
of the frequency, but depended only on the radial distance from the black hole horizon 
through ao. 

8.2 High Frequency Transverse Modes 
8.2.1 Electron-positron Plasma 

There exist four high frequency electromagnetic modes for the electron-positron plasma. 
These are illustrated in Figs. 5 and 6. Fig. 5 shows two modes which are similar and 
both show a large amount of growth. The bottom left mode is clearly a growth mode 
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Figure 7: Top: Real part of the growth mode and of the growth and damping mode for 
the high frequency transverse wave in the case of the electron-ion plasma. Bottom: Left: 
Imaginary part of the growth mode. Right: Imaginary part of the growth and damping 
mode. 

for all the frequencies, but the bottom right mode is growth for higher frequency and is 
damped very close to the black hole horizon for lower frequency. Thus near the horizon, it 
appears that energy is no longer fed into wave mode by the gravitational field but begins 
to be drained from the waves. The bottom left mode, shown in Fig. 6, is damped mode 
but the bottom right mode is growth for lower frequency as a) — > 1 and Oq —>■ 0. These 
four modes are similar to the three modes of Buzzi et al. [57] with larger damping and 
growth rates. 

8.2.2 Electron-ion Plasma 

Four high frequency modes are found to exist as shown in Figs. 7 and 8. The bottom 
left mode in Fig. 7 is a growth mode but the bottom right mode shows both growth and 
damping. For each of these modes, the solution become unstable for ao — and oj ^ 1. 
This may simply because the solution is too close to a resonant frequency. The damping 
and growth rates are larger than the corresponding modes for Schwarzschild case. Fig. 
8 shows two modes, one of which shown in bottom left is a purely damped, while the 
bottom right one is both damping and growth. These modes are stable for all frequencies 
and at all distances from the horizon. Unlike the corresponding Alfven modes, the growth 
and decay rates obviously depend on frequency. The four modes coalesce with the three 
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Figure 8: Top: Real part of the damping mode and of the damping and growth mode 
for the high frequency transverse wave in the case of the electron-ion plasma. Bottom: 
Left: Imaginary part of the one damped mode. Right: Imaginary part of the other mode 
showing both damping and growth. 

modes of Buzzi et al. [17] for the Schwarzschild case. 

9 Concluding Remarks 

The prime concern of this study has been exclusively the investigation, within the local 
approximation, of Alfven and high frequency transverse electromagnetic waves in a two- 
fluid plasma in the vicinity of the Reissner-Nordstrom-de Sitter black hole. 

We derive the dispersion relation for the Alfven and high frequency electromagnetic 
waves and solve it numerically for the wave number k. In the zero gravity limit, our study 
gives the results of the special relativistic case [53] where only one purely real mode was 
found to exist for both the Alfven and high frequency electromagnetic waves. In contrast 
to the special relativistic case, new modes (damped or growth) arise in the present work 
due to the black hole's gravitational field. We have found four Alfven and four high 
frequency modes either for electron-positron or electron-ion plasma. There is a variation 
of damping and growth rates of all the modes by some orders of magnitude with respect 
to the Schwarzschild case [17]. 

For the Alfven waves, the damping and growth rates are clearly independent of fre- 
quency for both the fluid components, but are solely dependent only on the radial distance 
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from the horizon as denoted by the mean- value of the lapse function a^. On the other 
hand, for the high frequency electromagnetic waves the damping and growth rates are 
dependent on both frequency and radial distance from the black hole horizon. So, we 
may conclude that the damping modes demonstrate the energy drain from the waves by 
the gravitational field and the growth rates indicate that the gravitational field is, in fact, 
feeding energy into the waves. 

The study of this paper gives the result for the lukewarm black hole when = 
M^. Moreover, the result can be specialized for the interesting cold, charged Nariai, and 
ultracold black holes by suitably choosing the black hole parameters as described in the 
introduction. In the limit £^ — oo the study provides the result for the RN [57] black hole, 
while for i"^ oo, = 0, the result reduces to that for the Schwarzschild black hole |47] . 
The result for the Schwarzschild-de Sitter black hole [58] is obtained if one sets = 0. 
Charged black holes in de Sitter space have interesting wormhole aspects [HU |59j. In 
gravitational theory wormholes are bridges between different parts of the universe. There 
exists a class of extreme cases which appears to be classically and quantum mechanically 
stable [521 EH]. The RNdS black hole is therefore interesting in a broader context. Our 
study of this paper is thus well motivated. 
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